Fracture of a viscous liquid 
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When a viscous liquid hits a pool of liquid of same nature, the impact region is hollowed by the 
shock. Its bottom becomes extremely sharp if increasing the impact velocity, and we report that the 
curvature at that place increases exponentially with the flow velocity, in agreement with a theory 
by Jeong and Moffatt. Such a law defines a characteristic velocity for the collapse of the tip, which 
explains both the cusp-like shape of this region, and the instability of the cusp if increasing (slightly) 
the impact velocity. Then, a film of the upper phase is entrained inside the pool. We characterize 
the critical velocity of entrainment of this phase and compare our results with recent predictions by 
Eggers. 

PACS numbers: 47.15.Rq, 47.15. Gf, 47.55. Dz 



We all observed that diving in a swimming-pool pro- 
duces air bubbles, while entering slowly in the same pool 
does not: the motion of a solid penetrating a liquid at a 
high speed largely deforms the surface of this liquid, and 
induces air entrainment. This phenomenon often defines 
the maximum speed of coating of a solid, which enters a 
bath before leaving it coated. Generally air entrainment 
must be avoided because of the resulting bubbles (which 
are all the more harmful since the coating solution gener- 
ally contains surfactants, and thus is likely to transform 
into a foam). Different experiments using either fibers or 
plates as solids indeed prove the existence of a threshold 
(in velocity) for air entrainment , but there is today no 
theoretical picture for understanding quantitatively this 
phenomenon. 

Quite similarly, when a viscous liquid enters a pool of 
same nature (as it occurs for example for the stem wave 
of a moving ship) , air entrainment may occur (producing 
bubbles, which ultimately burst - a serious cause of noise 
for a ship) . Such a bubble production was observed for a 
liquid jet impinging a bath of same nature, as reported by 
Lin and Donelly, who studied experimentally the mini- 
mum velocity of air entrainment as a function of different 
characteristics of the jet (diameter, viscosity or surface 
tension) 0. This question is important in many indus- 
trial processes where a viscous liquid (typically molten 
glass, metal or polymer) is poured inside a mould ~ and 
there again, the formation of bubbles is detrimental. A 
remarkable geometry for studying this phenomenon, and 
more generally the phenomenology of the regimes of de- 
formation of the liquid bath, was proposed by Joseph 
0, Q . Two counter-rotative cylinders are partially im- 
mersed in a bath, whose upper surface is observed. At a 
low speed of rotation, the interface between the cylinders 
is hardly deformed, but increasing this speed leads to the 
formation of a cusp-like shape. This device thus produces 
a singularity at a free interface, and more generally the 
observation and description of similar singularities has 
been a very active field for the past years, either theo- 
retically [5j or experimentally - let us quote for example 



selective withdrawal 0], the pinch-off of a jet 0, drops 
running-down on a tilted sohd @ or jet eruptions 0]. 

Surface tension naturally opposes the formation of 
such singularities at a liquid/fluid interface, and for the 
problem of the cusp between two counter-rotative cylin- 
ders, Jeong and Moffatt calculated in a classical paper 
the (non-zero) radius r of the tip, in the limit of small 
Reynolds numbers j^. They derived an analytical ex- 
pression, which can be summarized using a physical ar- 
gument proposed by Hinch In this two-dimensional 
geometry, the tip can be seen as a thin (hemi-)cylinder, of 
radius r. The viscous drag (per unit length) on a cylinder 
scales as rjV/ ln(c?/r), denoting d as an external length, -q 
as the liquid viscosity and V as the flow velocity around 
the cylinder. On the other hand, the surface tension 7 
draws this cylinder (as it does on the edge of a free sheet 
of liquid) with a force 27 (there again, per unit length). 
Balancing both these forces yields as a tip radius: 

r ~ (iexp(— Ca) (1) 

where Ca is the capillary number (Ca = rjV/^). 

This law first helps to understand why a tip forms 
around a well-defined velocity: the cusp becomes very 
sharp when V is of the order of 7/??, a logical quantity 
since both these parameters play antagonist roles in the 
formation of the tip. But it also predicts that r essen- 
tially vanishes above this velocity, and Eggers pointed 
out that this should lead to the crack of the tip, because 
of the lubrication pressure of the upper fluid, forced to 
flow in this very confined region jlQl] . He showed that 
beyond a critical radius for the tip, the pressure in the 
cusp breaks its stationary shape. Thus, this critical ra- 
dius Tc must increase with the viscosity 770 of the upper 
fluid, and Eggers found that Tc should scale as {ijo/vy^^ 
[13 • Together with Moffatt's result (Eq. [TJl, this yields 
as a threshold velocity Vc of entrainment of the upper 
phase: 

-ln(r//ryo) (2) 
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We could have expected (naively) that entrainment of 
the upper fluid would occur when the viscous forces over- 
come capillary ones {i.e. Vc ^ j/v)- Eggers' theory (Eq. 
O shows that taking into account the flow of the upper 
phase in the cusp should induce a logarithmic correction 
to this simple statement. As a consequence, the thresh- 
old velocity should depend on the viscosity of the upper 
phase, i.e. on its nature. This might sound a posteriori 
quite logical, but remarkably, this dependence is expected 
to be extremely weak. 

Here, we investigate experimentally the different sit- 
uations related to this problem. We flrst describe the 
macroscopic deformation of the interface, due to the in- 
coming flow of liquid. Then, we focus on the smaller 
scale of the tip. Finally, we characterize the condition of 
entrainment of the upper phase. 

We used a horizontal cylinder of radius R = 2 cm, half- 
immersed in a bath of glycerol (or in some cases silicone 
oil), of density p and viscosity rj around 1 Pa.s (Figure^l. 
The whole device was set in a transparent tank, and the 
upper fluid was either air or light oil (of density po < p 
and viscosity rjo <^ rj). The interfacial tension between 
the two fluids, denoted as 7, was measured for each ex- 
periment. The cylinder could be rotated at a controlled 
velocity f2, and the motion filmed with a digital camera. 
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FIG. 1: Experimental set-up. A solid cylinder is placed at the 
interface between a viscous dense liquid and a lighter fluid. 
The deformation of the interface between both fluids is mon- 
itored as a function of the speed of rotation of the cylinder. 

A viscous liquid film of thickness h (of typically 5 mm, 
thus thick enough to be independant of the nature of the 
roller surface) is dragged out of the bath. This layer im- 
pinges the bath at a velocity V ~ fl{R + h), and deforms 
the free surface at this place by a depth L. We denote 
r as the radius of curvature of the bottom of this hollow 
region. Figure[21shows the different shapes obtained with 
a silicon oil (77 = 0.97 Pa.s), and increasing the speed of 
the cylinder. 

It is observed that both the characteristic lengths L 
and r are dynamic in essence: the depth of immersion L 
increases with the speed of cylinder, while the tip radius 
r decreases. Moreover, these lengths are very different 
from each other: while the depth of immersion is around 
1 cm, the tip has a radius of curvature much smaller 
than 1 mm. This radius can become very small when the 
speed is large, as first reported by Joseph et al. Q in a 
similar situation. But we do not observe as they did a 
clear threshold velocity for the formation of a cusp - both 




FIG. 2: Four snapshots showing the bath deformation for 
an increasing speed of rotation of the cylinder {V is 6 cm/s, 
10 cm/s, 17 cm/s and 30 cm/s, going from left to right). The 
lower liquid is a silicon oil of viscosity 0.97 Pa.s and the upper 
fluid is air. The bar in the last picture indicates 1 cm. 

the lengths L and r varying continuously as a function 
of the cylinder velocity. 

We first focus on the macroscopic length L, namely the 
depth of immersion of the tip into the bath, and report 
in Fig. Olthe variation of L as a function of the impacting 
velocity V, for a bath of silicon oil (77 ~ 0.97 Pa.s) . 
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FIG. 3: Depth of immersion of the tip as a function of the 
velocity of the incoming liquid. The bath is a viscous silicon 
oil (77 = 0.97 Pa-s) with air above. The straight line indicates 
a slope 0.5. 

The depth L of immersion of the tip is found to increase 
with V. For speeds larger than 10 cm/s, L is several mil- 
limetres, which is significantly larger than the capillary 
length = y^j/Apg « 1.5 mm (denoting Ap as the 
difference of density between both fluids). We conclude 
that capillary forces can be neglected in this regime of 
large velocities. The viscous friction which pulls down- 
ward the tip can be written dimensionally riV/L^, while 
gravitational forces per unit volume (which pull it up- 
wards) are Apg. Balancing both these forces yields: 

L ^ ivV/Apg)'/' (3) 

Equation 121 defines some kind of viscous length (on the 
model of n^^ defined just above, replacing capillary ef- 
fects by viscous ones). It describes fairly well the data in 
Figure O where L is observed to increase as the square 
root of V. Some deviations appear at small V : the de- 
formation then is small {L ~ k~^), and capillary effects 
must be considered. Since they oppose the formation 



3 



of the hollow region, they indeed make L smaller than 
predicted by Eq. 13 

We then focused on the small scale of the deforma- 
tion, namely the tip itself, and took macrophotographs 
of this region. A numerical camera was placed behind an 
optical microscope lighted by a classical white lamp, and 
the whole device first calibrated using micrometric fibers. 
The tip was observed to be slightly deeper close to the 
container walls (which allowed us to make conveniently 
our measurements), but remains far from it (about 15 
mm because of the lubrification of the wall. As reported 
in Figure 01 the tip has a typical radius smaller than 
100 ym, which decreases strongly with the flow velocity. 




FIG. 4: Macrophotograph of the tip region (the bar indicates 
200 Jim) , for a tip of air in a bath of glycerol. The flow veloc- 
ities are respectively 14 cm/s and 22 cm/s, and the sharpness 
of the tip is observed to be highly dependant on this velocity. 

Analyzing the contour of the interface close to the tip 
allows a determination of the curvature: we can fit this 
contour by a parabolic function y = ax^ , from which we 
deduce the curvature r — l/2a. We display in Figure 
|5l the measured tip radii as a function of the capillary 
number Ca ~ l"i for air and oil tips formed in a bath 
of glycerol. 
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FIG. 5: Tip radius as a function of the capillary number 
Ca — rjV/"f for air tips formed in a bath of glycerol (circles, 
rj = 0.91 Pa-s, 7 = 55 mN/m), and oil tips in a bath of glycerol 
(squares, 77 = 0.53 Pa-s, 770 = 0.5 mPa-s, 7 = 22 mN/m). The 
data are deduced from macrophotographs such as the ones in 
igure H The plot is semi-logarithmic, and the equation of 
the straight fltting line is r = 993 exp( — 1.91Ca) (expressed 
in micrometers). 

The tip radius is found to decrease in a continuous 
way, as a function of the flow velocity V. More precisely, 
its logarithm varies linearly with y - in agreement with 
Moffatt's law (Eq. on nearly two orders of magnitude 



for the tip radius. This exponential behavior helps un- 
derstanding why the cusp-like shape apparently sets for 
a well-defined velocity, of the order of 7/77, although the 
actual variation of the radius is continuous with V. But 
it also suggests that slightly above this velocity, the cusp 
should become extremely sharp. This is not the case: a 
small increase of the roller velocity leads to the destruc- 
tion of the tip, and the upper phase then is entrained, as 
reported in Figure Note that the entrainment occurs 
simultaneously all along the tip, showing the 2D nature 
of the phenomenon. 




FIG. 6: Above a well-defined velocity Vc, a thin yet visible 
sheet of the upper fluid (here air), appearing in black in the 
image and stressed by an arrow, comes out of the tip and 
enters the viscous bath (here glycerol). The velocity V is 
32 cm/s, while the threshold Vc for air entrainment is 30 cm/s. 

We denote Vc as the threshold velocity and Cuc as 
the corresponding capillary number. Our experiments 
consisted in determining Cqc- Firstly, we varied the bath 
viscosity 77, using mixtures of water and glycerol and air 
as the upper fluid. 77 was measured for each solution. The 
corresponding data are displayed in black in Figure H 
Secondly, we used pure glycerol for the bath, and various 
light oils (such as hexane, hexadecane and silicone oils) 
as the upper phase. We measured 7 for each oil (it was 
found to be quite constant, around 26 ± 1 mN/m). The 
corresponding critical capillary numbers are represented 
by white points, in Figure |7| where all our results are 
plotted as a function of the viscosity ratio rjo/i], in a 
semi-logarithmic scale. 

We clearly observe an effect of the viscosity 770 of the 
upper fluid (open points): the larger 770, the smaller the 
critical capillary number. Nevertheless, this dependency 
is weak, since it appears to be well described by a loga- 
rithmic law. It is also limited: if the upper fluid is too 
viscous (typically 770/77 > 0.1), the lower liquid is not en- 
trained on the roller. For a given upper fluid, varying the 
bath viscosity also leads to a (weak) variation of the 
threshold in capillary number - which thus is not simply 
given by a value of order unity. All the data lie on the 
same semi-logarithmic plot, in agreement with Eggers' 
model (Eq. [3), for a variation on flve orders of mag- 
nitude of the viscosity ratio rio/rj. However, the slope 
deduced from the fit is —0.22, larger than —1 the value 
predicted by Eq. El 

On the whole, these experiments may help clarifying 
the question of the stationary pattern which sets as a 
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FIG. 7: Critical capillary number Ca^ = V^c/'y above which 
the upper phase is entrained in the viscous bath as a function 
of the ratio between both viscosities (semi- logarithmic scale). 
The full points are obtained using water-glycerol mixtures 
(250 mPa-s < rj <1500 mPa-s)with air above, while the open 
points correspond to glycerol (77 — 900 mPa-s) with above 
different light oils (0.3 mPa-s < r]o <83 mPa-s). 

viscous liquid hits itself, in a two-dimensional geometry. 
If the impinging velocity is large enough (i.e. for capil- 
lary numbers typically larger than 0.1), two very different 
dynamic lengths can be defined: a global one, which is 
the depth of liquid hollowed by the impact, and which 
results from a balance between viscous effects and grav- 
ity; and a local one, at the bottom of this hollow region 
which can be expressed as the radius of curvature at this 
place. As predicted by Moffatt j5|], this radius decreases 
exponentially with the capillary number. This law de- 
fines a characteristic velocity for which the hollow region 
becomes very sharp, i.e. takes a cusp- like shape, as first 
reported by Joseph 0. At the same time, it predicts 
a collapse of the tip radius for finite velocities, which 
implies the crack of the tip (entrainment of the upper 
phase), as postulated by Moffatt Q and shown theoreti- 
cally by Eggers 0| . We indeed report here such a crack, 
and show that the critical velocity for which it occurs is 
mainly governed by the bath viscosity, but also depends 
(in a much weaker manner) on the viscosity of the upper 
phase - in agreement with Eggers' theory. 



Different studies could naturally complement this one. 
Firstly, it would be worth understanding (both experi- 
mentally and theoretically) the thickness of the film of 
air, in the regime of entrainment. The selection of this 
thickness could remind the Landau-Levich problem (se- 
lection of the thickness of the film adhering to a solid 
withdrawn out of a bath), but the problem here is purely 
dynamical, which raises an interesting (and open) ques- 
tion (as in selective withdrawal, for the radius of the emit- 
ted filament) . The problem of the stability of this film, 
which eventually breaks in bubbles, also remains to be 
studied. Another issue concerns the practically impor- 
tant question of water, for which the threshold velocity 
is much higher, which would require to study the infiu- 
ence of inertia in all the different processes described all 
along this study. Secondly, we observed a deviation to- 
wards the law displayed in Figure [3 using a very vis- 
cous silicone oil as a bath (instead of glycerol), and air 
above. Then, for 77 > lO^rjQ, we measured a critical cap- 
illary number much higher, typically of the order of 5 to 
10 (and increasing rapidly as decreasing the ratio rjo/rj). 
A possible cause of this discrepancy could be the poly- 
meric nature of the liquid (for which the chains are all 
the longer since the viscosity is high), which could lead 
to non-Newtonian effects close to the tip - and thus to 
changes in the threshold of air entrainment. Using new- 
tonian silicone oils of smaller viscosity (with either air or 
ethanol above) indeed yields data in very close agreement 
with the ones displayed in Figure Finally, it would be 
useful to connect this study with similar ones in other di- 
mensionality. In particular, Cohen stressed that for selec- 
tive withdrawal {i.e. one-dimensional tip forming when 
sucking locally a liquid above an interface) , the viscosity 
of the fluid formi ng th e tip does not affect the value of the 
cut-off curvature [III . This could be due to the fact that 
the Stokes law on a (hemi-)sphere does not contain any 
logarithmic dependency, as it does for a (hemi-)cylinder, 
which finally emphasizes that the effects reported here 
are two-dimensional in essence. 
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